Confinement: a real-time visualization 
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Due to the mechanism of confinement, as known from quantum chromodynamics, it is difficult to 
observe individual particles carrying fractional quantum number [e.g. quark with fractional electric 
charge). A condensed matter example of fractionalized particles is spinous in quasi-one-dimensional 
spin systems, which are domain walls in the background of Neel configurations carrying spin-|. 
Using the time-evolving block decimation algorithm, we visualize the nontrivial spinon dynamics 
induced by the confine mechanism in a two- leg spin-^ ladder. It can be illustrated by a simple single- 
particle picture of Bloch oscillation, not only qualitatively but also quantitatively. We propose the 
experimental realization and the real time detection of the spinon dynamics in the ultra-cold boson 
systems of ^^Rb. 
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I. INTRODUCTION 

In quantum chromodynamics, quarks are fundamental 
particles carrying fractional electric charges. However, 
interactions between quarks grows linearly with distance 
due to the mechanism of confinement arising from the 
SU{3) color gauge theory. Consequentially, quarks are 
confined into color singlet bound states of baryons and 
mesons, and thus are difficult to observe. In various con- 
densed matter systems, there also exist excitations carry- 
ing fractional quantum numbers, such as the spinon and 
holon excitations in ID Luttinger liquids^, the solitons 
with half-fermion charge in one dimensional conducting 
polymers^, and the quasi-particle and hole excitations in 
fractional quantum Hall systems^, and the monomer ex- 
citations in the dimer model in the triangular lattice^. 
In these systems, fractional excitations are deconfined. 
Some of them have been experimentally detected. 

The phenomenon of confinement emerges in quasi- ID 
strongly correlated systems such as spin ladders. The 
"quark" in this case is known as spinon, which is the do- 
main wall interpolating between different ordered regions 
and usually carries fractionalized spin (spin-1/2). Take 
an antiferromagnetic spin chain with Ising-like exchange 
anisotropy for an example, the ground state in this case 
would be a Neel state with two-fold degeneracy, while 
the spinon can be considered as the elementary excita- 
tion separating these two degenerate ground states with 
opposite staggered magnetization. In a single spin chain, 
no confinement occurs because a pair of spinous could 
be separated as far from each other as possible without 
costing energy. For a two-leg ladder system, even an 
infinitesimal interchain coupling would induce a poten- 
tial which increases linearly with distance between the 
spinons^^. Therefore a pair of fractionalized spinon ex- 
citations would be confined into an integral spin-1 excita- 
tion: a magnon. Recently, a finite temperature confined- 
deconfined crossover has been observed in neutron scat- 
tering experiments for a weakly coupled ladder material 
CaCu203^. The signature of the confinement of spinous 



can be observed from the energy absorption spectrum for 
spin- flips at various wavevectors. In this paper, we reex- 
amine this old concept from a novel perspective, which 
enables us to visualize the confinement directly through 
the non-equilibrium dynamics in a cold atom system. 

Due to the low dissipation rate and the long coherence 
times, ultracold atoms in optical lattices have opened 
exciting possibilities for studying non-equilibrium quan- 
tum dynamics of many-body systems^iSrJ^. Qn the other 
hand, it also provides a perfect platform to reexamine 
classic concepts in condensed matter or particle physics 
from a different perspective. One example is the spin- 
charge separation, which plays a central role in strongly 
correlated systems^^^^^. Recently, Kollath et al^^ have 
used the adaptive t-DMRG^^ to study the time evolu- 
tion of a ID fermionic Hubbard model in real time, and 
observed the splitting of local perturbation into separate 
wave packets carrying charge and spin. In this paper, we 
study the time evolution of a pair of spinous in a two-leg 
spin ladder model, and find that the confinement mech- 
anism would lead to nontrivial dynamics of the spinous, 
which enables us to visualize confinement in real time. 



II. CONFINEMENT IN A TWO-LEG SPIN 
LADDER 

Our departure point is a two- leg spin-^ ladder with 
easy axis anisotropy along the z-axis defined as 
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(1) 



where a = 1,2 is the leg index; || and ± denote the 
couplings along the leg and across the rung, respectively. 
The couplings within legs are antiferromagnetic with the 
Ising anisotropy {J^ \\ > 2 J^^ y > 0). The rung coupling is 
chosen as Ising-like and ferromagnetic. The experimental 
realization of Eq.([T]) in ultra-cold atom systems will be 
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FIG. 1: Confinement on the two- leg spin-| ladder with Ising 
anisotropy described by Eq.([T|). (a) The classic ground state 
with Neel ordering along the leg and ferromagnetic ordering 
along the rung, (b) The classic configuration with spinoni 
and spinon2 located at the same position, (c) The mismatch 
of the locations of spinoni ^2 gives rise to the energy cost lin- 
early increasing with the mismatch length. 



discussed later. Below we will focus on the case of Ji y 7^ 
J2J1 which gives rise to the nontrivial dynamics induced 
by spinon confinement. 

Because of the Ising anisotropy and the ferromagnetic 
rung coupling, the ground state is a Neel state within 
each leg and with the spins of each rung aligned as shown 
in Fig. [U (a)^^. Introducing two spinous in each leg re- 
spectively located at the same position cost the domain 
wall energy £^5 = ^{Ji y + ^2 ||) shown in Fig. [1] (b). 
The separation of spinoni and spinon2 will lead to a mis- 
match within the rungs between spinoni and spinon2. 
This induces a finite string tension between spinous that 
grows with distance, leading to the confinement^, as 
shown in Fig. [I](c). The transverse spin couplings Ja^\\ 
generate spinon dynamics. In the absence of the rung 
coupling, two spinous move independently with veloci- 
ties Va cx ^a,||^^- III the case of Ji y 7^ J2J1, spinoni and 
spinon2 propagate with different velocities, thus becomes 
apart eventually. The effective potential between spinous 
is known as^: 
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where r is the average distance between these two spinon; 
a is the lattice constant; {G\Sf ^Sf 2\G) is taken over the 
ground state. 



III. CONFINEMENT-INDUCED DYNAMICS OF 
SPINONS 

In our numerical simulation, we use the open bound- 
ary condition (OBC) and prepare the initial state as both 



FIG. 2: Time evolution of the spinon dynamics with = 
with the spinon positions marked by the red cross. The 
spinous in the 1st (left) and 2nd (right) leg move indepen- 
dently. 



spinous located at the left end of the ladder, i.e., between 
the first and second sites. Right after the beginning of 
the evolution, the spinous will be rebounded by the left 
boundary and propagate towards the left end with a ve- 
locity proportional to y. The length of the ladder is 
set to L = 20, long enough for the time scales simulated 
to exclude finite-size effects. Utilizing the time-evolving 
block decimation (TEBD) algorithm^^, we study the time 
evolution of the many-body wavefunction from the initial 
state. Total conservation is used to reduce the com- 
putational effort. In the course of real time evolution 
we take the truncation dimension % = 80 and time step 
Ar = 0.05. The convergence is checked by taking larger 

The calculated time evolution for the case of = 
is shown in Fig. 2 in which two spinous are decoupled. 
The parameter values are taken as Ji y = J, J2,|| = 0.5 J, 
and J-^ y = ^2 II =5^- To visualize the spinon dynamics, 
we present the time-evolution of the configuration of the 
expectation value of Sz- The spinous are located at the 
bonds connecting two sites with the same sign of Sz- 
During the evolution, two spinous become separated due 
to their different velocities determined by Ja^\\- For a 
clear presentation of the spinon, we define the rectified 
magnetization S^^a for each leg as 



Si.a = -Sla X (-l)^ 



(3) 



The location of spinoni or spinon2 is determined at the 
bond across which S^^a changes the sign. The time evo- 
lution of the spatial distributions of Si^a=i,2 at = 
is depicted in Fig. [3] (a). Both legs exhibit the propa- 
gation of spinous at uniform speeds. The speed of the 
first leg is larger than the that of the second one because 
^1,11 > ^2,11 . 

Now we turn on the coupling along the rung, a con- 
finement potential emerges between two spinous. It acts 
like a "constant" force connecting these two spinous. If 
we view the spinous as particles, the dynamics of the 
spinous can be qualitatively understood from a simple 
two-body picture as follows. Two particles with different 
initial velocities and masses move in a ID lattice system 
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FIG. 3: The time evolution of the spatial distribution of 
Si,a=i,2 for (a) Ji = and (b) Jl = J. They exhibit the 
propagations and oscillations of the spinous in the first (left) 
and second (right) legs. 



and interact with each other via a linear potential. Each 
particle oscillates around the center of mass (COM) of 
the system, which continues to do uniform linear motion 
until it reaches the boundary. This simple picture could 
be verified qualitatively by our numerical simulation of 
the spinon dynamics. The time evolution of spatial dis- 
tribution of Si,a=i,2 at = J is presented at Fig. [3](b). 
An oscillation of the spinoni is observed, while that of 
spinon2 is not as clear as spinoni. The effective mass of 
spinon2 is larger than that of spinoni^ thus its oscillation 
amplitude is too small to be visible. 

To further verify the two-particle picture quantita- 
tively, we choose J2J1 = 0, which means spinon2 is lo- 
calized in its initial position at the boundary. In this 
case, the two-body problem reduces to a one-body prob- 
lem of spinoni in a lattice system under a static mag- 
netic field, which is provided by the rectified spins in the 
second leg. The time evolution of a spinon within a ferro- 
magnetic spin chain under a constant magnetic field has 
been studied previously, the dynamics turns out to be a 
perfect Bloch oscillation^^. In our case, the confinement 
'constant' force provided by the static spinon2 plays a 
similar role of the constant electric field {E ex J^), and 
we expect a similar Bloch oscillation dynamics. It implies 
that the frequency of the oscillation of spinoni should 
be proportional to the interaction strength: 1/T ex J^, 
where T is the oscillation period. This relation can be 
verify numerically, as shown in Fig. HI where we can find 
a linear relation between 1/T and J]_ for both Ji y = J 
and 1.5 J. Simulation times were long enough to observe 
the stability of the oscillation phenomenon. 

In a perfect Bloch oscillation, the frequency is deter- 
mined by the strength of the external field (J^ in our 
case), and should be independent of the band width (Ji y 
in our case). However, we find that the oscillation fre- 
quency not only depends on J^, but is also slightly de- 
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FIG. 4: The relation of 1/T v.s. J]_ with parameter values of 
J2,\\ = 0, and Jiji = J (squares) or 1.5 J (circles). 



pendent on Ji y, especially for large values of Jj_. There 
are two possible factors contributing to the imperfection 
of Bloch oscillation. One is the boundary effect. Con- 
sidering the quantum nature of spinous (wave packets), 
the boundary begin to influence the dynamics even be- 
fore the centers of spinous reach the boundary. For larger 
J^, the amplitude of the Bloch oscillation A is smaller 
(A oc Jij|/J^), which indicates the motion of the spinoni 
is confined to the boundary, therefore the boundary effect 
becomes more obvious. This can explain why the devia- 
tion becomes larger for bigger Ji^\\/ Jj_' The other reason 
comes from deviations from the linear relations of Eq. 
(2), which can be considered as a zero-order approxima- 
tion of the exact interaction between spinous. Actually, 
on a long length scale, it is possible that higher order 
nonlinear modifications such as l/|r| emerge due to the 
quantum nature of spinous^. 



IV. EXPERIMENTAL REALIZATION 

For an experimental realization of our two-leg spin lad- 
der system in an ultracold atomic superlattice system, 
one first has to create Neel-ordered chains along the x- 
direction^^, for example by using ^^Rb atoms with two 
internal spin states |T = l^mp = il)- The necessary 
superexchange interactions of the ultra cold atoms can 
be achieved by the use of a superlattice syste m^^i^^i^^~ — . 
In the optical lattice, a typical energy scale for the su- 
perexchange coupling J is in the order of kHz, which 
implies a ms time-scale real-time dynamics in our case. 
By modifying the bias between neighboring lattice sites, 
both ferromagnetic and antiferromagnetic superexchange 
interaction can be realized^^; the orientation of the setup 
would be as shown in Fig. [5l The two-leg ladder struc- 
ture is achievable by applying another laser along the y- 
direction to construct a second superlattice structure. In 
addition, the easy-axis anisotropy is realizable by tuning 
the interaction U^^ = Un 7^ U^i'^^'^^ or by a periodi- 
cally modulated lattice^. Local addressing techniques 
as developed recently^ would be used to create domain 
walls, e.g. by a global spin-fiip on one half of the system, 
as well as for the read-out procedure. These techniques. 
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while advanced, are now experimentally established. The 
arguably most difficult part of the experiment would be 
to adjust the exchange interaction separately along two 
different legs. As mentioned above, the exchange in- 
teraction can be adjusted by superlattices^^, but to do 
this separately on legs 1 and 2, an additional superlat- 
tice structure is needed. In fact, antiferromagnetic spin 
chains have already successfully been simulated, albeit in 
a different context, namely the Bose-Hubbard model in a 
tilted optical lattices, where a magnetic domain wall has 
been observed^. 




FIG. 5: Sketch picture for the experimental setup of our two- 
leg spin ladder system. 



V. DISCUSSION AND CONCLUSION 

In the context of condensed matter physics, the study 
of domain wall motion in magnetic nanowires has at- 
tracted considerable theoretical and experimental atten- 



tions recently due to its potential industrial applications. 
Manipulation and control the domain wall dynamics in 
magnetic nanowires is known to play an important role 
in nanomagnetism^- — . Most of the theoretical stud- 
ies in this field are based on the Landau-Lifshitz-Gilbert 
(LLC) equation^^'^^. Our work is different from previous 
work in three aspects: (i) we provide a quasi-exact calcu- 
lation of the dynamics of the strongly correlated system 
by TEBD instead of a semi-classic mean-field approxima- 
tion; (ii) different from condensed matter system where 
the domain wall dynamics is closely related with dissipa- 
tion rate, there is little energy dissipation in our system; 
(iii) we account for the strong confined interaction be- 
tween domain walls, which leads to nontrivial dynamics. 

Whether in condensed matter or ultracold atom 
physics, the control and manipulation of the particles 
or quasiparticles in the quantum system are of central 
interest. In this paper, confinement is visualized in real 
space and its properties are studied via a non-equilibrium 
process, showing a possible way to investigate and ma- 
nipulate the motion and interaction of quasiparticles in 
a dissipationless strongly correlated system. Our work 
could be helpful for the study of phenomena involving 
the domain wall dynamics, such as Walker breakdown^ 
and current-driven domain wall motion^^. 
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supported by DFG through FOR801. 



^ D. Senechal, arXiv:cond-mat/9908262 (1999). 

2 A. J. Heeger, S. Kivelson, J. R. Schrieffer, and W. -P. Su 

Rev. Mod. Phys. 60, 781 (1988). 
^ R. B. Laughlin, Rev. Mod. Phys. 71, 863 (1999). 

R. Moessner, S. L. Sondhi, Phys. Rev. Lett 86, 1881 (2001). 
^ E. Dagotto and T. M. Rice, Science 271, 618 (1996). 
^ D. G. Shelton, A. A. Nersesyan, and A. M. Tsvelik, Phys. 

Rev. B 53, 8521 (1996). 

M. Greiter, Phys. Rev. B 66, 054505 (2002). 
^ M. J. Bhaseen and A. M. Tsvelik, |arXiv:cond-mat /0409602 
(2004). 

^ B. Lake, A. M. Tsvelik, S. Notbohm, D. A. Tennant, T. G. 
Perring, M. Reehuis, C. Sekar, G. Krabbes and B.Biichner, 
Nature Phys. 6, 50 (2010); M. Greiter, Nature Phys. 6, 5 
(2010). 

^° T. Kinoshita ,T. Wenger, D. S. Weiss, Nature. 440, 900 
(2006). 

K. Winkler et aL, Nature. 441, 853 (2006);A. J. Daley et 
a/., arXiv:cond-mat/0608721 
^2 S. Trotzky, Y. A. Chen, A. Flesch, I. P. McCulloch, U. 
Schollwock, J.Eisert, LBloch,arXiv:cond-mat/l 101. 2659 
K. A. Al-Hassanieh, F. A. Reboredo, A. E. Feiguin, 1. 
Gonzalez, and E. Dagotto, Phys. Rev. Lett. 100, 166403 
(2008). 



C. L. Hung, X. B. Zhang, N. Gemelke, and C. Chin, Phys. 
Rev. Lett. 104, 160403 (2010). 

D. Chen, M. White, C. Borries, B. DeMarco, arXiv: 
11034662 

C. Kasztelan, S. Trotzky, Y. A. Chen, 1. Bloch, I. P. Mc- 
Culloch, U. Schollwock, G. Orso, Phys. Rev. Lett. 106, 
155302 (2011). 

E. H. Lieb, F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968); 
G. Baskaram, Z. Zou, P.W. Anderson, Solid State Com- 
mun. 63, 973 (1987). 

C. Kollath, U. Schollwock, and W. Zwerger, Phys. Rev. 
Lett. 95, 176401 (2005); A. Kleine, C. Kollath, 1. P. Mc- 
Culloch, T. Giamarchi, and U. Schollwock, Phys. Rev. A 
77, 013607 (2008). 

T. Giamarchi, Quantum Physics in One Dimensions, 
(Clarendon Press, Oxford, 2004). 

A. J. Daley et aL, J. Stat. Mech. (2004) P04005; S. 
R. White and A. Feiguin, Phys. Rev. Lett. 93, 076401 
(2004); U. Schollwock, Rev. Mod. Phys. 77, 259 (2005); 
U. Schollwock, Ann. Phys. (N.Y.) 326, 96 (2011). 

D. Gobert, C. Kollath, U. Schollwock and G. Schiitz, Phys. 
Rev. E 71, 036102 (2005). 

G. Vidal, Phys. Rev. Lett. 91, 147902 (2003); Phys. Rev. 
Lett. 93, 040502 (2004). 



5 



Z. Cai, L. Wang, X. C. Xie, U. Schollwock, X.R. Wang, 
M. Di Ventra ,and Y.R Wang, Phys. Rev. B 83, 155119 
(2011). 

S.Trotzky, P. Cheinet, S. Foiling, M. Feld, U. Schnor- 

rberger, A.M. Key, A. Polkovnikov, E.A. Demler, M.D. 

Lukin, and I. Bloch, Science. 319, 295 (2008). 

Y. A. Chen, S. Nascimbene, M. Aidelsburger, M. Atala, S. 

Trotzky, I. Bloch, Phys. Rev. Lett. 107, 210405 (2011). 

S.Trotzky, Y. A. Chen, U. Schnorrberger, P. Cheinet, and 

I. Bloch, Phys. Rev. Lett. 105, 265303 (2010). 

J. J. Garcia-Ripoll, M. A. Martin-Delgado, and J. L Cirac, 

Phys. Rev. Lett. 93, 250405 (2004). 

L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 
91, 090402 (2003). 

A. B. Kuklov and B.V. Svistunov, Phys. Rev. Lett. 90, 
100401 (2003). 



C. Weitenberg, M. Endres, J. F. Sherson, M. Cheneau, P. 
SchauB, T. Fukuhara, L Bloch, S. Kuhr, Nature 471, 319 
(2011) 

J. Simon, W. S. Bakr, R.C. Ma, M. E. Tai, P. M. Preiss, 
M. Greiner, Nature 472, 307(2011) 

S. S. P. Parkin, M. Hayashi, and L. Thomas, Science 320, 
190 (2008). 

G. S. D. Beach, et al., Nat. Mater. 4, 741 (2005). 

M. Klaui, C. A. F. Vaz, J. A. C. Bland, W. Wernsdorfer, 

G. Faini, E. Cambril, L. J. Heyderman, F. Nolting, and U. 

Riidiger, Phys. Rev. Lett. 94, 106601 (2005); 

N. L. Schryer and L. R. Walker, J. of Appl. Phys. 45, 5406 

(1974). 

X. R. Wang, P. Yan , J. Lu, C. He, Annals of Physics 324, 
1815 (2009). 



